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Abstract 

The Bergman-Milton bounds provide limits on the effective permittivity of a com- 
posite material comprising two isotropic dielectric materials. These provide tight 
bounds for composites arising from many conventional materials. We reconsider the 
Bergman-Milton bounds in light of the recent emergence of metamaterials, in which 
unconventional parameter ranges for relative permittivities are encountered. Specif- 
ically, it is demonstrated that: (a) for nondissipative materials the bounds may be 
unlimited if the constituent materials have relative permittivities of opposite signs; (b) 
for weakly dissipative materials characterized by relative permittivities with real parts 
of opposite signs, the bounds may be exceedingly large. 
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1 Introduction 

Increasingly, new materials which exhibit novel and potentially useful electromagnetic re- 
sponses are being developed [1, 2}. At the forefront of this rapidly expanding field lie meta- 
materials [3]. These are artificial composite materials which exhibit properties that are 
either not exhibited by their constituents at all, or not exhibited to the same extent by their 
constituents. With the emergence of these new materials — which may exhibit radically 
different properties to those encountered traditionally in electromagnetics/optics — some 
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re-evaluation of established theories is necessary. A prime example is provided by the recent 
development of metamaterials which support planewave propagation with negative phase 
velocity [4]. The experimental demonstration of negative refraction in 2000 prompted an 
explosion of interest in issues pertaining to negative phase velocity and negative refraction 
[5,6]. 

The process of homogenization, whereby two (or more) homogeneous constituent ma- 
terials are blended together to produce a composite material which is effectively homoge- 
neous within the long-wavelength regime, is an important vehicle in the conceptualization 
of metamaterials [7]. The estimation of the effective constitutive parameters of homoge- 
nized composite materials (HCMs) is a well-established process [8], aspects of which have 
been revisited recently in light of the development of exotic materials that exhibit properties 
such as negative phase velocity. For example, it was demonstrated that two widely used 
homogenized formalisms, namely the Maxwell Garnett and Bruggeman formalisms, do not 
provide useful estimates of the HCM permittivity within certain parameter regimes [9]. The 
Maxwell Garnett estimates coincide with the well-known Hashin-Shtrikman bounds [10] on 
the HCM permittivity. While the former are commonly implemented for both dissipative 
and nondissipative HCMs, the later were derived for nondissipative HCMs. 

In view of the limitations of the Maxwell Garnett and Bruggeman formalisms within 
certain parameter regimes, we explore in this communication the implementation of the 
Bergman-Milton bounds [11, 12, 13] for these parameter regimes. To be specific, we consider 
the homogenization of two isotropic dielectric constituent materials with relative permittiv- 
ities e a and e&. We explore the regime in which the parameter 4 

is negative-valued, as this is where the Maxwell Garnett and Bruggeman estimates are not 
useful [9]. Notice that the definition (1) caters to the possibility that only one of e a or e& e K, 
as might arise for a metal-in-insulator HCM, for example. In fact, the 5 < regime which 
occurs for metal-in-insulator HCMs [14, 15] is highly pertinent to the homogenization of 
HCMs which support planewave propagation with negative phase velocity. 

Let us note that, although the 5 < regime has been discussed in the past in the context 
of the Bergman-Milton bounds [12, 15, 16], the discussion on the inadequacy of those bounds 
has been brief. Amplification is needed because of the possibility of fabricating negatively 
refracting composite materials [17, 18], for example. 

2 Bergman— Milton bounds 

Two bounds on the effective relative permittivity e e of the chosen composite material were 
established by Bergman [11, 19, 20, 21] and Milton [12, 22]. We write these as BM a and 

4 Re (e a ,b) and I m {ca,b) denote the real and imaginary parts of e a ,b, respectively; R and C denote the sets 
of real and complex numbers, respectively. 
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BM^. In terms of a real-valued parameter 7, these are given by (see eqn. (24) in [12]) 



and 



BM a ( 7 ) = f a e a + f b e b 



bm^( 7 ) = <; k + h. 



fafb{Cb 



3 [ye a + (1 - 7) e b ] ' 
2 fafb (e a - eb) 2 



e a e& 3e a e{, [e b 7 + e a (1 - 7)] 



(e a ,e b GC), (2) 



(e a ,e b EC), (3) 



where / ai (, denotes the volume fraction of the constituent material with relative permittivity 
e a> b, and f a + fb = 1- For the bound BM Q the parameter 7 takes the values (1 — / ) /3 < 
7 < 1 — /a/3, whereas for the bound BM^ the parameter 7 takes the values 2 (1 — f a ) /3 < 
7<l-2/ /3. 

The Bergman-Milton bounds (2) and (3) are related to the two Maxwell Garnett esti- 
mates of the HCM relative permittivity [8, 14] 



MG a = e b + 
MG,3 = e a + 



e a + 2e 6 - / a (e a - e 6 ) ' 
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(5) 



e b + 2e a - f b (e b - e a )' 

The Maxwell Garnett estimates represent the extension of the Hashin-Shtrikman bounds 
[10] into the complex-valued permittivity regime. For nondissipative HCMs, the Maxwell 
Garnett estimates coincide with the Bergman-Milton bounds when the parameter 7 attains 
its minimum and maximum values; i.e., 
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In view of our particular interest in homogenization scenarios for which 5 < 0, we note 



that 
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(7) 



for nondissipative mediums. Thus, there exist 

(i) a volume fraction f a G (0, 1) at which MG Q is unbounded for all values of 5 < —2, and 

(ii) a volume fraction f a G (0, 1) at which MG/3 is unbounded for all values of 5 G (—1/2, 0). 
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3 Numerical illustrations 



Let us now numerically explore the Bergman-Milton bounds, along with the Maxwell Gar- 
nett estimates, for some illustrative examples of nondissipative and dissipative HCMs. The 
parameter 5, defined in (1), is used to classify the two constituent materials of the chosen 
HCMs. We begin in §3.1 by considering nondissipative HCMs. While these do not represent 
realistic materials, they provide valuable insights into the limiting process in which weakly 
dissipative materials become nondissipative. Furthermore, they provide a useful yardstick 
in the evaluation of dissipative HCMs, which are considered in §3.2. 

3.1 Nondissipative HCMs 

We begin with the most straightforward situation: nondissipative HCMs arising from con- 
stituent materials with 5 > 0. In Figure 1, the Maxwell Garnett estimates MG a and MG^ 
(which in this case are identical to the Hashin-Shtrikman bounds) are plotted against 
f a G (0, 1) for e a = 6 and e = 2. The Bergman-Milton bound BM Q is given for f a e 
{0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9}. The corresponding plots of BM^) with 7 overlies 
that of BM a . The Bergman-Milton bounds are entirely contained within the envelope con- 
structed by the Maxwell Garnett estimates. 

Let us turn now to the nondissipative scenario wherein 5 < 0. In Figure 2, the the Maxwell 
Garnett estimates MG a and MG/3 are presented as functions of f a for e a = —6 and e& = 2. 
The Bergman-Milton bound BM Q is given for f a e {0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9}. 
The corresponding Bergman-Milton bound BM^ is plotted in Figure 3. In consonance with 
(6) and (7), we see that MG a becomes unbounded as f a — > 0.25. It is clear that MGp < 
BM Q < MG a for f a < 0.25, whereas MG a < BM^ < MG^ for f a > 0.25. For f a > 0.25, the 
Bergman-Milton bound BM Q lies outside both Maxwell Garnett estimates MG Q and MG^, 
and similarly BM^ lies outside both Maxwell Garnett estimates MG Q and MG^ for f a < 0.25, 
although the relations (6) still hold. 

3.2 Dissipative HCMs 

We turn to homogenization scenarios based on dissipative constituent materials; i.e., e a ,b £ C. 
Let us begin with the 5 > scenario. In Figure 4, the homogenization of constituents charac- 
terized by the relative permittivities e a = 6+0. 3i and e& = 2+0. 2i is illustrated. In this figure, 
the Maxwell Garnett estimates on complex- valued e e are plotted as f a varies from to 1. The 
Bergman-Milton bounds, which are graphed for f a e {0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9}, 
are fully contained within the Maxwell Garnett envelope. That is, we have MG^ < BMq, ^ < 
MG a for all values of f a . 

Now we consider dissipative constituent materials with 5 < 0. In Figure 5, the homoge- 
nization of constituent materials given by e a = — 6 + 3i and e& = 2 + 2i is represented. The 
Maxwell Garnett estimates are plotted for f a e (0, 1), whereas the Bergman-Milton bounds 
are given for f a e {0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9}. As is the case in Figure 4, BM^ lies 
entirely within the envelope constructed by MG a and MG^. We see that BM a > MG^ for all 
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values of f a ; but, for mid-range values of f a , BM a slightly exceeds MG Q for certain values of 
the parameter 7. 

As the degree of dissipation exhibited by the constituent materials is decreased, the 
extent to which BM Q exceeds MG a is increased. This is illustrated in Figure 6 wherein the 
homogenization is repeated with e a = — 6 + i and e b = 2 + 2i/3. As in Figure 4, the Maxwell 
Garnett estimates are plotted for f a e (0, 1), while the Bergman-Milton bounds are given for 
f a E {0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9}. The Bergman-Milton bound BM^ lies within the 
Maxwell Garnett envelope for all values of f a , but substantial parts of BM a lie well outside 
the envelope of the two Maxwell Garnett estimates. 

The behaviour observed in Figures 5 and 6 is further exaggerated in Figure 7, where 
the homgenization of constituent materials with e a = — 6 + 0.3i and e b = 2 + 0.2i is rep- 
resented. The Maxwell Garnett estimates are plotted for f a e (0, 1); for reasons of clarity, 
the Bergman-Milton bounds are plotted only for f a e {0.1,0.3,0.5}. The Maxwell Garnett 
estimates are exceedingly large and the Bergman-Milton bounds are larger still. 

Finally, let us focus on the scenario referred to in the introduction, namely the homog- 
enization of a conducting constituent material and a nonconducting constituent material, 
wherein 5 < 0. Suppose we consider constituents characterized by the relative permittivities 
e a = — 6+3i and e b = 2. In Figure 8 the Maxwell Garnett estimates are plotted for f a e (0, 1), 
whereas the Bergman-Milton bounds are given for /„ e {0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9}. 
As we observed in Figure 6, the Maxwell Garnett envelope does not contain substantial parts 
of the Bergman-Milton bound BM a , whereas the BM^ bound lies entirely within the envelope 
constructed from the two Maxwell Garnet estimates. 

4 Discussion and conclusions 

The Bergman-Milton bounds, as well as the Maxwell Garnett estimates, are valuable for 
estimating the effective constitutive parameters of HCMs in many commonly encountered 
circumstances. However, the advent of exotic new materials and metamaterials has led to 
the examination of such bounds within unconventional parameter regimes. It was recently 
demonstrated that the Bruggeman homogenization formalism and the Maxwell Garnett ho- 
mogenization formalism do not provide useful estimates of the HCM permittivity when the 
relative permittivities of the constituent materials e a and are such that [9] 

(i) Re (e ) and Re (e b ) have opposite signs; and 

(ii) |Re(e 0>6 ) | > |Im(e a , b ) |. 

Similarly, we have demonstrated in the preceding sections of this communication, that the 
Bergman-Milton bounds do not provide tight limits on the value of e e within the same 
parameter regime. 

We note that if the real parts of e a and e b have opposite signs, but are of the same order 
of magnitude as their imaginary parts, then the Bergman-Milton bounds are indeed useful, 
and they then lie within the envelope constructed by the Maxwell Garnett estimates. 
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Figure 1: The MG Q (thick dashed line) and MG^ (thin dashed line) estimates of e e plotted against f a 
for e a — 6, eb = 2. The vertical solid lines represent the variation of the Bergman-Milton bound BM a 
with 7 for f a e {0.1(a), 0.2(6), 0.3(c), 0.4(d), 0.5(e), 0.6(/),0. 7(5), 0.8(/i), 0.9(i)}; and these coincide with the 
corresponding variation of BM^ with 7. 
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Figure 2: The MG Q (thick dashed line) and MG^ (thin dashed line) estimates of e e plotted against f a 
for e a = —6 and e = 2. The Bergman-Milton bound BM Q is plotted as the vertical broken lines for 
f a e {0.1(a), 0.2(6), 0.3(c), 0.4(d), 0.5(e), 0.6(/), 0.7(g), 0.8(h), 0.9(t)}. 
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Figure 3: As Figure 2 but with BM/3 (vertical solid lines) in place of BM 
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Figure 4: The MG Q (thick dashed line) and MG/3 (thin dashed line) estimates in relation to Ree e 
and Ime e as f a varies from to 1, for e a = 6 + 0.3i, = 2 + 0.2i. The Bergman-Milton bounds 
BM Q (thin broken dashed lines) and BMp (thin solid lines) in the top diagram are plotted for f a G 
{0.1(a), 0.2(6), 0.3(c), 0.4(d), 0.5(e), 0.6(/), 0.7(g), 0.8(/i),0.9(i)}. The bottom diagram shows the Bergman- 
Milton bounds in greater detail but for f a = 0.5(e) and f a = 0.6(f). 
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Figure 5: The MG Q (thin dashed line) and MGp (thick dashed line) estimates in relation to 
Ree e and Ime e as f a varies from to 1, for e a = —6 + 3i and — 2 + 2i. The Bergman- 
Milton bounds BM Q (thin broken dashed lines) and BM^ (thin solid lines) are plotted for f a G 
{0.1(a), 0.2(6), 0.3(c), 0.4(d), 0.5(e), 0.6(/), 0.7(g), 0.8(h), 0.9(i)}. 
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Figure 7: As Fig. 4 but for e a = -6 + 0.3« and £& = 2 + 0.2i The Bergman-Milton bounds are plotted for 
/« G {0.1(a), 0.3(c), 0.5(e)}. 
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Figure 8: As Fig. 4 but for e a = -6 + 3i and e b = 2. 
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